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Abstract. We outline a general algorithm for computing an explicit 
model over a number field of any curve of genus 2 whose (unpolarized) 
Jacobian is isomorphic to the product of two elliptic curves with CM by 
the same order in an imaginary quadratic field. We give the details and 
some examples for the case where the order has prime discriminant and 
class number one. 

1 Motivation 

Let E\ , E2 be two elliptic curves defined over Q with complex multiplication by 
an order O of an imaginary quadratic field K. We are interested in finding explicit 
models for curves C defined over Q whose (unpolarized) Jacobian is isomorphic 
to Ei x E2. In this paper we propose a general algorithm for this purpose and 
give details only for the following special case where we have have carried them 
out: Ei = E2, O is the ring of integers of K = Q(\/—N), N = 3 mod 4 prime 
and O has class number one. Our special case consists of finitely many curves, 
up to isomorphism; the algorithm produces models over K for them. 

It is a general fact due to Narasimhan and Nori [NN] that there are only 
finitely many principal polarizations on a given abelian variety up to isomor- 
phism. Hence, for a fixed O there are only finitely many isomorphism classes of 
the curves we want; their number was calculated by Hayashida and Nishi [HN]. 

For a similar question in the case of abelian surfaces with complex multipli- 
cation by a quartic field see [vW] . 

Our interest in this problem arose in connection with a generalization to 
genus 2 of the singular moduli formulae of Gross and Zagier [GZ] for the norm 
of the difference of j-values of CM elliptic curves. (This generalization will be 
the subject of a separate publication.) As an illustration, consider the genus 2 
curve C determined by 

y 2 = f(x) = 6- 3 h(x)V{x), 
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where 



h(x) = (7144V 3 163 - 151790)a; 3 + (129789^/^X63 + 1752597)x 2 
+ (-47481V^163 + 510153)a;+ (-1596^/^163- 37250) , 

and 

h l (x) =x 3 h(-l/x) 

(bar denoting complex conjugation of the coefficients). The unpolarized Jacobian 
of C is isomorphic over Q to the product of two elliptic curves with CM by the 
ring of integers of K = Q(V— 163). Let 

D = 2- 12 disc (/) = (2 • 3 2 • 5 • 7 • 11 • 17 • 19 • 23) 12 . 

Then we have 




where 

/24 4 6 \ 
Q{m) = m* 4 55 1 I rn 
\ 6 1 83/ 

is a certain positive definite ternary quadratic form of level 163 associated to 
C and the (finite) sum is over m € Z 3 such that (163 — Q(m))/4 is a positive 
integer. In particular every rational prime I dividing D is smaller than 163/4 
and inert in K. 

The significance of the number D is that C has bad reduction only at primes 
dividing D. Note that over Q the Jacobian of C has good reduction everywhere 
but C does not; at primes diving D, C reduces to two elliptic curves crossing at 
a point. 

Another source of interest in the problem is the fact, which I learned from K. 
Lauter, that the reduction of the curves C provides genus 2 curves over certain 
finite fields with maximal number of rational points (see §5 for an example). 
In this regard, the more interesting problem is the analogous one for curves of 
genus 3 for which we hope to exhibit in the near future an algorithm similar to 
the one sketched here. 

2 Outline of the algorithm 

We start by giving an outline of the main steps of the general algorithm and 
then give details for our special case in the next sections. 

Step 1. Find period matrices for the polarized Jacobians. 

Step 2. Given a non-split period matrix obtained in step 1 compute a model 
for the corresponding curve. 



The first step is purely algebraic and only requires computations with rational 
numbers; it involves the calculation of representatives for ideal classes of certain 
orders in a quaternion algebra (see [HN] for more details). What we need to do 
is describe explicitly the finitely many principal polarizations on E\ x E 2 up to 
equivalence. 

The second step relies on the following explicit version of Torelli's theorem 
for curves of genus 2 due to Bolza and Klein. Let 

H 2 = {Z E C 2x2 | Z* = Z, Im(Z) positive definite} 

be the Siegel upper- half space of rank 2. Let Z E TL 2 be a period matrix of 
a principally polarized abelian surface which is not the product of two elliptic 
curves with the product polarization. A theorem of Torelli guarantees that Z 
arises from a curve of genus 2, unique up to isomorphism. Here is a way of 
recovering the curve from the period matrix Z . 

Let fz(ui,u 2 ) E C[ui,u 2 ] be the leading term in the Taylor expansion of 

J J 0^, v (u, Z) , u=(ui,u 2 ) 
ip.v) odd 

about the origin, where for fi, v € {0, 1} 

0„,„(u, Z) := Yl e WZm e 2W ( u +H , u = ( Ul , u 2 ) E C 2 , Z E H 2 

meZ ! +i(i 

is the theta function with characteristics (see [Mu]). Then the canonically po- 
larized Jacobian of the hyperelliptic curve over C determined by the equation 

y 2 = .fz{x,l) 

corresponds to Z. (There are six theta functions with odd characteristics and 
hence fz is a sextic, i.e. homogeneous of degree 6.) 

The difficulty in applying this formula is to know how to normalize the sex- 
tic fz properly to guarantee that its coefficients are algebraic integers as well 
as finding similar expressions for its Galois conjugates. In general, this would be 
accomplished by an application of Shimura's general reciprocity law. We would 
obtain rapidly convergent series giving the minimal polynomials of these coef- 
ficients. Since the coefficients of the minimal polynomials of the coefficients of 
the sextic are in Z, truncating the series would then allow us to compute them 
exactly. We show how this works for our special case in §4. 

3 Principal polarizations 

From now on we assume that O is the ring of integers of K = Q(\/— ~N) C C, 
N = 3 mod 4 prime and the class number of O is 1. Hence, E\ = E 2 = E, with 
E isomorphic to C/O over C. 



The principal polarizations of Ex E up to isomorphism correspond to positive 
definite unimodular Hermitian forms of rank 2 over O up to GX2(0)-equivalence. 
In order to find a set of representatives of these Hermitian forms we will exploit 
the happy accident that since we assume N = 3 mod 4 the quaternion algebra 
B = ( ~ 1 q" JV ) (up to isomorphism the unique quaternion algebra over Q ramified 
only at N and oo) contains Q(i). This allows us to convert the question to that 
of finding Hermitian forms over Z[i] of discriminant — N up to equivalence and 
this is quite simple. Here is how it works. 

Consider in B the order 

R = Z + Zi + t}-{\ + j) + Zil(l + j), i 2 = -l,j 2 = -N. 

R is a maximal order in B with a natural embedding of O sending \/—N to 
j. The rank 2 unimodular Hermitian forms arising from polarizations of E x E 
correspond to rank 1 left i?-modules. 

Since R also has an embedding of Z[i] (sending i to i) we may associate to 
a left R- module a rank 2 Hermitian form <P over Z[i]. We can give <P as a triple 
(a, b, c) with a, c € Z >0 and 6 e Z[z], where <£(u, w) = 2cum + 6mT + 6uw + 2cvv. 
It is not hard to see that this form has discriminant bb — 4ac = — N. 

It will be more convenient to work with SL2 rather than GL2 equivalence and 
to avoid duplications we consider only forms <P — (a,b,c) with 6=1 mod 2. The 
above discussion establishes a 1-1 correspondence between principal polarizations 
on E x E, up to SX2 (O)-equivalence, and positive definite binary Hermitian 
forms <P = (a,b,c) over Z[i] of discriminant — N, up to SL 2 (Z [^-equivalence. 

Let H be the hyperbolic 3-space 

H = {w = (x,y,t) € R 3 I t > 0} , 

which we will think as embedded in the Hamilton quaternion algebra H by 
(a;, y, t) ^ x + iy + jt (here i, j are the usual basis of H with i 2 = j 2 = — 1 and 

ij = —ji). If ^ e SL 2 (C) then it is not hard to check that 

w 1 ► (aw + b)(cw + d)^ 1 

sends H to H defining an action of 5X2(0) on H and in particular, an action of 
SX 2 (Z[«]). This last action has a very simple fundamental domain, whose closure 
is given by w = (x,y,t) e H with x 2 + y 2 + t 2 > 1, x < 1/2, y < 1/2, < x + y. 

We can associate to a form # the point w = (b + y/~Nj)/2a e H and we call 
# reduced if w lies in the fundamental domain. The action of SL 2 (Z[i]) on H 
mimics that on Hermitian forms. Every form is SL 2 (Z[i])-equivalent to a unique 
reduced form. 

The situation is in fact very analogous to that of positive definite binary 
quadratic forms over Z and, as in that case, it is easy to write an algorithm that 
lists all reduced forms <P of a given discriminant (we do not really need TV to be 
prime or class number 1 for this) . Here is a brief sketch (all of this is classical 
going back to Hermite [He, I, p. 251]). 



Input : N = 3 mod 4 

For < r, s < r odd, s even 

Set m := (r 2 + s 2 + N)/4 
For a|m, max(r, s) < a < ^fm 

Add (a,r + zs,m/a) to List 

Add (a, —r + si,m/a) to List unless 

a = m/a or r = a or s = a or s = 

Output : List 

As an example, we give in table 1 the list of reduced forms of discriminant —163. 



Table 1. Reduced Hermitian forms over Z[i] of discriminant —163 



In general, the number of <£'s is the class number n of B [Ei] , which can be 
given in terms of N as follows (a formula valid for any prime N = 3 mod 4) 



is a period matrix corresponding to the associated principal polarization on 
E x E. 

4 Bolza-Klein sextics 

The product polarization on E x E corresponds to the reduced form <P = 
(1, 1, (N + l)/4) in the principal class; hence, forms not in the principal class 
correspond to curves. 



(1,1,41) 
(2,1 + 2i,21) 
(3,±l + 2i,14) 
(6, ±1 + 2i,7) 
(4,±3 + 2i,ll) 
(6,±5 + 2i,8) 
(5,±l + 4i,9) 
(7, ±5 + 6i,8) 



Finally, given a <P ~ (a, 6, c 





Given a form <P = (a, b, c) not in the principal class we define the associated 
normalized Bolza-Klcin sextic as follows. 



where r\ is Dedekind's eta function. It satisfies the following properties. 

— The SLi{G) class of depends only on the SL 2 (Z[i]) class of <£. 

— /<p has coefficients in K and a 6 has coefficients in O. 

— The Igusa invariants [Ig] of are in Z and depend only on the SX 2 (Z [in- 
equivalence class of <P. 

The genus 2 curve 

C<p : y 2 = /<p(or, I) 

is then defined over if and, over the algebraic closure if of if in C, its Jacobian 
is isomorphic to E x E. 

Given a form <P = (a,b,c) let <P L := (a, —b,c). Suppose both <P and ^' are 
reduced. Then # and <P L are not 5i2(Z[i])-equivalent but they are (always) 
GL2(Z[z])-equivalent. The corresponding curves C<p 7 C<e>l are hence isomorphic 
over K; note that they are also complex conjugates of each other. Otherwise, 
curves C$ corresponding to different reduced forms are non-isomorphic. The 
involution l has a natural counterpart on the left i?-ideals in B and it turns out 
that the number of orbits of i is what is classically known as the type number 
of B [Ei]. Hence, there are t — 1 isomorphism classes of curves with Jacobian 
isomorphic to E x E, where t is the type number of the quaternion algebra B 
[HN]. 

Here is a table with the values of n and t for the primes N we are considering. 



Table 2. Type and class number of the quaternion algebra B 



N 


n 


t 


3 


1 


1 


7 


1 


1 


11 


2 


2 


19 


2 


2 


43 


4 


3 


67 


6 


4 


163 


14 


8 



Note that for N = 3 or 7 we only have the product polarization and hence 
there is no curve C with unpolarized Jacobian isomorphic to E x E in that case. 

Given a curve C defined over Q its field of moduli is the field F C Q char- 
acterized by the property: For every r G GaZ(Q/Q), C T is isomorphic to C if 



and only if t is the identity on F. Clearly isomorphic curves have the same field 
of moduli. Notice that F is the smallest field over which a curve isomorphic to 
C could be defined, but it is not in general a field over which it can be defined. 
In fact, for example, Shimura showed that no generic hyperelliptic curve of even 
genus has a model over its field of moduli [Sh Thm 3] . See [Me] for a discussion 
of this issue for curves of genus 2. 

For the curves C<p the field of moduli is Q (the field generated by the Igusa 
invariants [Ig]), but, in fact, most are not definable over Q; only those forms <P 
which are S'L2(Z[i])-equivalent to <P L give rise to curves definable over Q. 

To see this we note that by their very construction the period matrices lies 
in a certain real 3-dimensional cycle in TL2 considered by Shimura [Sh] . Namely, 
the cycle defined by 

(i^-K-ii)- 

If Az is the complex abelian surface corresponding to such a Z then there is an 
isomorphism 

A : Az — ► Az, with A o A = —id. 

(Applied to Z<|> this yields the fact that the curves Cp and are both isomor- 
phic and complex conjugate to each other as mentioned above). 

It follows that if Az has no automorphisms other than ±id then it has no 
model defined over its field of moduli. It is not hard to see that this holds for 
Z$, for every in the interior of the fundamental domain. 

5 Examples 

We end with an illustration of the above discussion, giving the outcome of algo- 
rithm when N — 43. The calculations were done using PARI-GP. The routines 
as well as the data for all cases is available at: 
http : //www.ma.utexas . edu/users/villegas 

The reduced forms <P of discriminant —43 are (1,1,11), (2,1 + 2i,6) and 
(3, ±1 + 2*. 4). 

1) For <5 = (2, 1 + 2i, 6) we obtain 

f (x, 1) = ~(- x 6 + i(-3 + 567 v / ^43) x A + i(3 + 567v^43) x 2 + 1) 
Its Igusa invariants are 

1728012 
93313728006 
186622271996 
2176943579975806271997 
2176782336000000000000 



J 2 = 

Je= - 
Js= - 
J 10 — 



(these were calculated using classical algorithms for invariants of a sextic follow- 
ing Mestre [Me]). 

As in the example of the introduction D = J w = 2~ 12 disc (/$) factors nicely 

D= (2 2 -3-5) 12 . 
This curve descends to Q; here is a model 
y 2 = x 6 + 24384 x 5 + 61311 x 4 + 585856 a; 3 + 813483 x 2 + 3214656 x + 1472877 . 

2) For $ = (3, 1 + 2i, 4) we obtain 

U(x,l) = j^((14>/=43- 160) x 6 +(42V^43+ 162) x 5 + 

(2247^/^43 - 159) x 4 +17021 x 3 + 
(22477^43 + 159) x 2 + (-42\/^43 + 162) x 
+14V^43+ 160) 

Its Igusa invariants are 

J 2 = 14333772 

J 4 = 7393823156166 

J 6 = 3726840435157546564 

J 8 = -312234946681873274015037 

J 10 = 7355827511386641000000000000 

and 

D = Ji = (2-3-5-7) 12 . 

As explained in §4, since ^ corresponds to a point in the interior of the 
fundamental domain the curve C<j> has no model over Q. Alternatively, we can 
see this following Mestre [Me]. When the curve has no extra automorphisms (i.e. 
its only automorphisms are the identity and the hyperelliptic involution), the 
obstruction to the curve being definable over its field of moduli (Q in our case) 
is given by a conic in P 2 

xMx 1 =0, x = (xi : x 2 : x 3 ) e P 2 , 

where M is a 3 x 3 symmetric matrix whose entries are certain invariants of the 
sextic; more precisely, the curve is definable over its field of moduli if and only 
if the conic has a rational point there. Explicitly, we have M = (rriij) with (we 
have actually simplified slightly the matrix given by Mestre) 

mn = 3J| - I6OJ4./2 - 3600J 6 

m 2 i = - J 4 Ji + 330J 6 J 2 + 160JJ 

m 3 i = - J 6 J| - 840J 6 Ji - 8OOOJ10 

m 22 = -25 J 6 Jf - 8 J| J 2 - 120 J 6 J 4 - 2000 J w ' 

m 32 = 67J 6 J 4 + 6OOJ10 J 2 + 90J| 

m 33 = -33J| J 2 - 100J 6 J| - 8OOJ10 J 4 



where J2, J4, Je, Js, Jw ar e the Igusa invariants. 
In our case we have 



mn = -21538723388574481387776 

mi2 = 24856361223852137345176064256 

mi 3 = -23971255400369899892885589544571136 

m 22 = -28732882146400381994651008552571136 

m 23 = 27776672840855638207256856144392139100416 

m 33 = -26987491534155851141341724256178812956900004096 

We easily verify that this conic has rational points everywhere locally except at 
the primes 43 and 00; in particular, it has no rational points. 

We should point out that the vanishing of the determinant of M precisely cor- 
responds to the curve having extra automorphisms. As with D, this determinant 
factors nicely 

dct M = -2 64 • 3 38 • 5 34 • 7 28 • 19 4 • 29 2 • 37 2 • 43 . 

Finally, let p be a prime which splits in K = Q(V— 43) as p = VP. The 
reduction of the curve C$ modulo V gives a smooth curve C of genus 2 over 
F p . We have verified that for all primes in the range 167 < p < 10000 such that 
Ap = a 2 + 43 for some a € N the curve C or its quadratic twist attains the 
maximum number of points possible, namely p+1 + 2 [y/2p ] (an improvement 
on Weil bounds due to Serre). 
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